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Abstract 

A new molecular simulation method for the calculation of vapour-liquid equilibria of mix- 
tures is presented. In this method, the independent thermodynamic variables are tem- 
perature and liquid composition. In the first step, one isobaric isothermal simulation for 
the liquid phase is performed, in which the chemical potentials of all components and 
their derivatives with respect to the pressure, i.e., the partial molar volumes, are calcu- 
lated. From these results, first order Taylor series expansions for the chemical potentials 
as functions of the pressure //^(p) at constant liquid composition are determined. That 
information is needed, as the specified pressure in the liquid will generally not be equal 
to the equilibrium pressure, which has to be found in the course of a vapour simulation. 
In the second step, one pseudo grand canonical simulation for the vapour phase is per- 
formed, where the chemical potentials are set according to the instantaneous pressure 
using the previously determined function ^ii{p"). In this way, results for the vapour pres- 
sure and vapour composition are achieved which are consistent to the given temperature 
and liquid composition. The new method is applied to the pure Lennard-Jones fluid, a 
binary, and a ternary mixture of Lennard-Jones spheres and shows very good agreement 
with corresponding data from the literature. 
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1 Introduction 



Due to the rapid increase of available computing power, molecular simulation is develop- 
ing to become a standard tool in chemical engineering and related areas. It can provide 
a better understanding of industrial processes and the underlying physics on the molecu- 
lar level. Furthermore, molecular simulation has the capability to predict phase equilibria 
more reliably than classical thermodynamic models. The calculation of vapour-liquid equi- 
libria by molecular simulation is a longstanding and important task. In the last 15 years 
a variety of methods for this purpose have been presented. 

There are, among others, thermodynamic scaling [T|2] . histogram reweighting [3P] . Gibbs- 
Duhem integration ^5ii6l| . NpT plus test particle method [TO||TT], various extensions of it to 
other ensembles [7|S1I^ . and, most widely used, the Gibbs Ensemble Monte Carlo method 
(GEMC) [T^PH] . A discussion of the different approaches can be found in the reviews 

mm- 

The GEMC specifies the thermodynamic variables temperature, global composition, and 
either global volume or pressure. Both phases are simulated simultaneously in separate 
volumes, between which real particles are exchanged in order to ensure phase equilib- 
rium. On the one hand, this set of specified thermodynamic variables is in many cases not 
convenient, i.e., when bubble points or dew points are needed, where temperature or pres- 
sure and the composition of one phase is specified. On the other hand, the simultaneous 
simulation of both phases has the disadvantage that fluctuations occuring in one phase 
influence the other phase directly. This can result in unfavourable dynamic behaviour of 
the simulation, e.g., close to the critical point fluctuations can lead to phase reversal in 
the two volumes. Furthermore, real particle exchange between the two phases is a major 
cause of statistical uncertainties. If the liquid is dense and/or consists of strongly inter- 
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acting particles, the insertion and deletion of real particles is awkward and perturbs the 
hquid causing high statistical uncertainties. 

In this work, the Grand Equilibrium method for the calculation of vapour-liquid equilibria 
of multi-component mixtures is presented. This method circumvents the above mentioned 
problems. The specified thermodynamic variables are temperature and liquid composition, 
two independent simulations for both phases are performed and no exchange of real par- 
ticles in the liquid phase needs to be done. 

The Grand Equilibrium method has some common features with the NpT plus test particle 
method for mixtures [TTfTB] . In fact, the treatment of the liquid phase is identical. The 
difference lies in the procedure for the vapour phase, where only one simulation has 
to be performed and no prior knowledge of the vapour-liquid equilibrium is neccessary. 
Furthermore, the present method is related to the method of lingerer et al. [T71IT8] . but 
avoids introducing direct coupling between the phases. Similar to the work of Esobedo 
[TUPU] . it employs pseudo ensemble simulations. 

The Grand Equilibrium method is tested in the present work on the pure Lennard- Jones 
fluid, a binary, and a ternary Lennard- Jones mixture and the results are compared to 
hterature data obtained by other methods. 
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2 Method 



The basic idea of the Grand Equihbrium method is to evaluate the chemical potentials 
in the liquid as functions of the pressure at a given temperature and liquid composition 
and to use these functions to set the chemical potentials in a pseudo grand canonical 
ensemble simulation for the vapour phase. The phase equilibrium conditions are directly 
introduced in the vapour simulation, which yields vapour pressure, vapour composition 
and other properties on the dew line. 

Let us consider a mixture consisting of an arbitrary number of components i = A, B, C, ... 
and a given set of intermolecular potentials between all species. Having set the independent 
thermodynamic variables temperature T and liquid composition x = xa, xb, xq, one 
simulation in the isobaric isothermal (NpT) ensemble is carried out at an arbitrarily 
chosen pressure p\ preferrably not too far from the vapour pressure p^j. This simulation 
yields the density = p(T,:s.,p'-) and the enthalpy = /i(T, x,p'). Applying fluctuation 
formulae in the NpT ensemble allows the derivatives of the density and enthalpy with 
respect to the pressure at constant temperature and constant composition to be obtained 
from that simulation. The derivative of the density is calculated from ^21] 



where V denotes the instantaneous volume, the total number of particles, k the Boltz- 
mann constant, /St the isothermal compressibility, and the brackets < > indicate averaging 
in the NpT ensemble. The derivative of the enthalpy is calculated from 



'dh\ _ 1_ 
9^ It,. " ^ 



.kT 

where H denotes the instantaneous enthalpy. 



^ ' < H ><V > - < HV > ]+ <V > 



(2) 
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For the present method it is essential to additionally calculate the chemical potentials 
/ij and the partial molar volumes f , of all components from this liquid phase simulation. 
This can be done, e.g., by Widom's test particle insertion p2j, according to which 



= ^f{T) + kT-\nxi-kT-\n<V- expi-iPi/kT) > /N, (3) 

where ipi denotes the potential energy of a test particle of species i and fil'^{T) the part 
of the ideal chemical potential which only depends on the temperature. It is convenient 
to define the residual chemical potential of component i as 

/i, = - fif{T)]/kT. (4) 



Widom's test particle insertion also yields the partial molar volumes 

<V'.expi-^,/kT)> 
= <V. eM-^./kT) > - < ^ > • 



It is important to note that the Grand Equilibrium method is not tied to Widom's test par- 
ticle insertion, which is known to yield high statistical uncertainties in very dense and/or 
strongly interacting fluids. In such cases, other more sophisticated methods to determine 
the chemical potentials, such as gradual insertion in extended ensembles P^25f26j . can 
be applied alternatively. 

The complete set of thermodynamic quantities that have to be calculated from one NpT 
simulation at given T, x, p'- in the liquid phase is p\ h'-, {dh/dpYq^^, /i', and f ■. 

With these quantities, the desired chemical potentials, as functions of the pressure, are 
obtained from first order Taylor series expansions [11] 



Mp) ~ /i- + v[/kT-{p-p'), (6) 
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with 



df4 
dp 



) 



T,x 



(7) 



It should be mentioned that in typical applications, the accuracy of equation ([6]) is not 
influenced significantly by neglecting the higher order terms, but rather by the accuracy 
of fil and f ■ achieved in the simulation. Furthermore, in a liquid phase not too close to 
the critical point the chemical potentials are only weak functions of the pressure. Q 

Having determined the fii{p) from the Taylor series expansions, according to equation ([6]) 
by one NpT simulation in the liquid mixture, one vapour simulation in a pseudo grand 
canonical ensemble (pseudo-/iVT) has to be carried out. In a common grand canonical 
ensemble (/uVT) simulation [27|28f29] . the temperature, the volume, and the chemical 
potentials of all components are set as fixed values. In the present pseudo-/iVT simulation 
the temperature and the volume are set in the common way. But instead of setting fixed 
values for the chemical potentials, they are set as functions of the instantaneous pressure 
in the vapour p"" using equation ([6]) to find iiiip"). Applying this procedure, it is ensured 
that equilibrium between the liquid phase and vapour phase is imposed. 

During a common /iVT simulation the chemical potentials are set through insertion and 
deletion of particles by the comparison between the resulting potential energy change 
and the desired fixed chemical potential. In the proposed pseudo-/iVT ensemble, in which 
the values of the chemical potential are found from the function /ij(p^), the acceptance 
probability Pace for particle insertion of species i writes as 



^ If equation Q is needed in the near-critical region and no good guess for is available, 
succesive simulations in the liquid phase can be necessary. 



Pacc{N^N + 1) = min 1, 



V 





Ni + l 
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where iVj denotes the number of particles of species i and U{N + 1) the configurational 
energy of the system with + 1 particles. The acceptance probability for particle deletion 
of species i writes as 

PacciN^N - 1) = min (l, y ■ exp{-/i,(j9") + [U{N) - U{N - 1)] /fcT}) . (9) 

Starting from an arbitrarily chosen low density state point, the vapour simulation is forced 
by this procedure to change its state to the corresponding phase equilibrium state point 
on the dew line. The vapour pressure p„ and the vapour composition on the dew line 
y = Ha, Vb, Vc, ■■■ are then simple ensemble averages in the vapour simulation in the 
pseudo-yuVT ensemble. The same holds for other properties on the dew line, such as dew 
density and dew enthalpy. 

Experience shows that the pseudo-yuVT simulation moves rapidly into the vicinity of the 
phase equilibrium state point. This process happens in the equilibration period. During 
the production period the pressure just fluctuates around its constant average. These 
inherent fluctuations in the pressure impose fluctuations on the set chemical potentials 
through equation (JH]). Such an interdependence characterizes pseudo statistical ensembles. 

Figure 1 gives an impression of a typical equilibration process obtained using the procedure 
described above for a binary Lennard- Jones mixture. After an initial equilibration of 1 000 
Monte Carlo loops in the canonical ensemble (NVT), a second equilibration in the pseudo- 
/iVT ensemble is started. It can be observed that the instantaneous values of pressure, 
density, and vapour composition reach the vicinity of equilibrium after only 1 000 Monte 
Carlo loops. Due to the experience with various mixtures including nonspherical and polar 
fluids [HT], it can be stated that this process is very robust and independent on the initial 
density and composition in the vapour simulation. 

With the knowledge of the vapour pressure, other thermodynamic properties of interest 
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on the bubble line can be calculated easily, e.g., the bubble density 

p' ^ p' + p'/3^(p.-p'). 



(10) 



the bubble enthalpy 



h' 



+ 



( 



dp 



) 



(11) 



and the chemical potentials in phase equilibrium 



picT 



(12) 



The present method can be applied for the calculation of dew points as well. Therefore 
the phases simply have to be inverted, i.e. one NpT simulation in the vapour phase and 
one pseudo-pV^T simulation in the liquid phase have to be performed. In this case it faces 
the awkward insertion and deletion of real particles in the liquid phase, but still has the 
advantage of independent simulations runs. 

Furthermore, the Grand Equilibrium method can be applied to the more challenging 
liquid-liquid equilibria as well. First results [BO] confirm the applicability and demonstrate 
high accuracy. 
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3 Applications 



In order to test the Grand Equilibrium method, it has been applied to the pure Lennard- 
Jones fluid, a binary, and a ternary Lennard- Jones model mixture. The Lennard- Jones 
intermolecular potential acting between two particles of species i and j is defined by 



Uij{r) 
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CTi 



(13) 



where r is the intermolecular distance and aij and Sij are the specific size and energy 
parameters. In the case of a Lennard- Jones mixture, the full set of size and energy pa- 
rameters for all interactions has to be specified. Unlike interaction parameters are often 
expressed on the basis of the like interaction parameters by 



cr. 



(a,, + a,)/2. 



(14) 



and 



(15) 



where ^j,- is a binary interaction parameter. 



3.1 Pure Lennard- J ones fluid 

The pure Lennard- Jones fiuid, being a very simple but realistic molecular model, has 
been studied extensively by many authors. Its vapour-liquid equilibria are known with 
high accuracy [3^33[34] , so that it is a favourable test fiuid here. 

In a pure fiuid the size and energy parameters are constant for all interactions, so that 
CTjj = a and Sij = e. All thermodynamic quantities are reduced in the usual way: temper- 
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ature T* = Tk/e, pressure p* = pcr^/e, density p* = pa^, configurational energy u* = u/e, 
configurational enthalpy h* = h/e, isothermal compressibility = Pt^/o'^, and residual 
chemical potential as given by equation (jlj). 

For a test of the Grand Equilibrium method on the pure Lennard- Jones fluid, no liquid 
simulations have been performed in the present work. As for the liquid phase simulations, 
only standard methods are needed, i.e. NpT ensemble and Widom's test particle insertion, 
so results from the literature [32] which are known to be accurate were used: Table 1 
comprises this data. 

Vapour simulations following the Grand Equilibrium method have been performed based 
on these results. The technical details can be found in the appendix. The vapour-liquid 
equilibrium results obtained with the Grand Equilibrium method are given in table 2. That 
table also contains the widely accepted results of Lotfi et al. [32j, which were obtained 
using the NpT plus test particle method [lO]. The statistical uncertainties of the present 
vapour-liquid equilibrium results on the dew line were determined by block averaging [35] 
of the vapour simulations. In this way, the influence of the uncertainties of the chemical 
potentials and the partial molar volumes from the liquid simulation is neglected. It should 
be emphasized that this leads to an underestimation. For a repeatability study, see section 
3.2. The statistical uncertainties of the vapour-liquid equilibrium results on the bubble 
line were determined by block averaging [33] of the liquid simulation and the utilisation 
of the error propagation law. 

It can be seen from table 2 that all data points for vapour pressure, bubble density, and 
dew density from the present work and those obtained by Lotfi et al. [32] agree within their 
very small statistical uncertainties. The uncertainties of the Grand Equilibrium method 
are lower throughout than those of the NpT plus test particle method. In our opinion, this 
is due to the above mentioned neglecting of the uncertainties of the chemical potentials 
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and the partial molar volumes from the liquid simulation, since this is an important source 
of error and the quality of the results of both methods should be similar. 

3.2 Binary Lennard- Jones mixture 

The Grand Equilibrium method was also tested on a binary Lennard- Jones model mixture, 
of which the vapour-liquid equilibrium is well known. The parameters of the investigated 
symmetrical binary model mixture A+B are cr^ = cr^, Sa = ^b^ and ^ab = 0.75, cf. 
equations ( fT3l) -(fT5l). The thermodynamic quantities are reduced in the same way as for 
the pure Lennard- Jones fluid, using the parameters of the A-A interaction, a a and Ea- 
For comparison to the data obtained in the present work, results of the following authors 
are used: Panagiotopoulos et al. [13] and Tsangaris et al. who applied the Gibbs 
Ensemble Monte Carlo method (GEMC) as well as Boda et al. [9] who applied the Grand 
Canonical Monte Carlo method (GCMC). 

Binary vapour-liquid equilibrium data is calculated for the two temperatures T* = 1 and 
T* = 1.15. Table 3 comprises the simulation results for the liquid phase. It should be 
pointed out that all liquid simulations on one isotherm have been performed at the same 
arbitrarily chosen pressure p'*. On the basis of these results, the vapour-liquid equilibria 
were determined by pseudo-yuVT simulations and are given in table 4. The technical details 
of the simulations can be found in the appendix. The statistical uncertainties indicated 
in table 4 were estimated with the same procedure as described above for the pure fluid. 

In order to investigate the statistical uncertainties of the Grand Equilibrium method 
appropriately, the vapour-liquid equilibrium for T* = 1.15 and xa = 0.05 have been 
calculated independently four times, cf. table 5. A better estimation of the uncertainties of 
the vapour-liquid equilibrium can be made by comparing these results. The uncertainties 
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are for vapour composition 6yA = 0.002, vapour pressure 6p* = 0.001, bubble density 
6p'* = 0.005, dew density 6p"* = 0.002, bubble enthalpy 6h'* = 0.03, and dew enthalpy 
6h"* = 0.02. These values are higher throughout than those obtained by the statistical 
analysis of the vapour runs, which is easily explained by neglecting of the uncertainties 
of the chemical potentials and the partial molar volumes from the liquid simulation on 
the vapour-liquid equilibrium. They are still considerably lower than those reported by 
Panagiotopoulos et al. [I3] for their method. 

Different comparisons between the results of the present work obtained using the Grand 
Equilibrium method and GEMC results of Panagiotopoulos et al. [I3] and Tsangaris et 
al. [SH] as well as GCMC results by Boda et al. [9j are given in figures 2-4. From the 
pressure-composition vapour-liquid equilibrium diagram in figure 2, it can be seen that 
the present data shows very little scatter, in the range of less than 0.01 in composition, 
and that it corresponds very well to the vapour-liquid equilibrium of the pure fluid B at 
XA = yA = 0. The agreement with the GEMC results P^|3S] is excellent, as almost all 
vapour-liquid equilibrium points agree with the GEMC results [THII^ within the statistical 
uncertainties. The agreement with the GCMC results [9] at the lower temperature is good, 
especially on the dew line. But at the higher temperature considerable deviations occur. 
For a given pressure, the GCMC [9] compositions on both the bubble and dew line seem 
to be too low by 0.02-0.03. 

In figure 3 a pressure-density vapour-liquid equilibrium diagram is given. Again, the 
present data shows very little scatter and is consistent to the results for pure fluid B, 
for which the saturated densities can be found as the lowest point of each isothermal 
branch. The agreement to the GEMC results p^f36] is excellent, as all points are within 
the statistical uncertainties of the GEMC results [T^IHU] . The dew densities of the GCMC 
method agree very well with the other data, but the bubble densities of [H] especially 
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at the higher temperature are somewhat too high. 

Figure 4 shows a pressure-energy density vapour-hquid equihbrium diagram, where the 
energy density is defined by the product of the configurational energy and the density. Due 
to the lack of appropriate hterature data for T* = 1.15, only the isotherm T* = 1 is shown. 
Once more, it is found that the energy density obtained with the Grand Equilibrium 
method shows very little scatter and corresponds very well to the limiting pure fiuid B, 
which is included as the lowest point of the present data. The agreement with GEMC 
|13|I5B] is again excellent. Consistent with the densities, the GCMC method [H] agrees well 
on the dew line, whereas on the bubble line, particularly for higher pressures, deviations 
are observed. 

3.3 Ternary Lennard- Jones mixture 

For a further test, the Grand Equilibrium method was applied to a ternary Lennard- Jones 
model mixture A+B+C, for which vapour-liquid equilibrium results from Tsang et al. [37J 
using the GEMC method, were available. 

That ternary mixture is defined on the basis of equations (fT^ - (fT^ by cr^ = cr^ = o'C; 
Eb = 0.75eA, = 0.15eA, and C,ab = (,ac = isc = 1- All thermodynamic quantities are 
reduced in the usual way, using the parameters of the A-A interaction, cr^ and Ea- 

In the first step, liquid simulations for different compositions have been performed at 
T* = 1 and p* = 0.2. The full set of these results is given in table 6. The vapour-liquid 
equilibrium data is calculated by pseudo-/iVT simulations in the vapour phase as described 
above, cf. table 7. The technical details of the simulations are given in the appendix. The 
lower bounds of the statistical uncertainties of the vapour-liquid equilibrium, as given in 
table 7, were estimated by the same procedure as described above for the pure fiuid. It 
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should be mentioned, that for each vapour-hquid equihbrium point only one simulation 
in the liquid phase and one simulation in the vapour phase is neccessary, regardless of the 
number of involved components. 

Figure 5 shows a ternary plot of the vapour-liquid equilibrium at T* = 1 and p* = 0.2. 
The results from the present work are compared to the GEMC results of Tsang et al. [37J. 
It can be seen that both data sets are in excellent agreement. 
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4 Conclusion 



A new Grand Equilibrium method for the calculation of vapour-liquid equilibria is pre- 
sented. It is aimed at mixtures, but can be applied to pure fluids straightforwardly. The 

independent variables are temperature and liquid composition which is favourable for 
bubble point calculations. 

Two subsequent independent simulation runs have to be performed, first in the liquid, then 
in the vapour, where standard simulation techniques are utilized. The method requires no 
prior knowledge of the vapour-lquid equilibrium, is easy to use, and performs reliably. 

The Grand Equilibrium method has been applied successfully to the pure Lennard-Jones 
fluid, a binary, and a ternary mixture. Comparisons to results from other methods show 
that it yields reliable results with very little scatter. 
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Appendix 

For the liquid phase, NpT molecular dynamics simulations were performed, using the 
algorithm of Andersen [SH] and applying periodic boundary conditions as well as the 
minimum image convention. A cut-off radius of Aa^ was introduced, considering the long 
range corrections. The equations of motion for 500 particles were solved with a fifth-order 
predictor-corrector method using a time step of 0.005 in usual units. Starting from a f.c.c. 
lattice, the equilibration period was 5 000 time steps, whereof the first 1 000 were carried 
out in the NVT ensemble. The length of the production run was 100 000 time steps. In 
order to calculate the chemical potentials and the partial molar volumes, Widom's test 
particle method [22] was utilised, by inserting 1 000 test particles after each time step. 
The statistical uncertainties of the simulation results were determined according to the 
method of Fincham [33] . 

For the vapour phase, pseudo-/xVT Monte Carlo simulations were performed, applying 
periodic boundary conditions and the minimum image convention. The cut-off radius was 
again Aa^ and the long range corrections were considered. The maximum displacement 
was set to 5% of the simulation box length, which was chosen to yield on average 300 to 400 
particles in the volume. After 1 000 initial NVT loops starting from a f.c.c. lattice, 9 000 
equilibration loops in the pseudo-yuVT ensemble were performed. One loop is defined here 
to be a number of attempts to displace particles equal to the actual number of particles 
plus three insertion and three deletion attempts. The length of the production run was 
100 000 loops. 
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Table 1 

Thermodynamic properties in the hquid phase of the pure Lennard-Jones fluid taken from Lotfi 
et al. [32] • The numbers in paranthesis denote the uncertainty in the last digit. 





P 






P* 






0.75 


0.00 


-5.69 


(3) 


0.8214 


(4) 


0.091 (3) 


1.00 


0.03 


-3.823 


(7) 


0.7018 


(4) 


0.28 (2) 


1.15 


0.06 


-3.200 


(4) 


0.6056 


(6) 


0.6 (1) 


1.25 


0.10 


-2.882 


(4) 


0.516 


(2) 


2.3 (5) 



21 



Table 2 

Comparison of vapour-liquid equilibria results of the pure Lennard-Jones fluid from the present 
work (Grand Equilibrium method) and from Loth et al. [32] (NpT plus test particle method). 
The numbers in paranthesis denote the uncertainty in the last digit. 







* 

Pa 






p'* 






p"* 






0.75 





00261 


(1) 





8216 


(4) 





00359 (1) 


present work 







00264 


(7) 





8216 


(4) 





0036 


(1) 


m 


1.00 





02500 


(4) 





7008 


(4) 





0296 


(1) 


present work 







0250 


(2) 





7008 


(4) 





0296 


(3) 


m 


1.15 





0597 


(2) 





6054 


(6) 





0728 


(4) 


present work 







0597 


(4) 





6055 


(7) 





0727 


(8) 


[32J 


1.25 





0971 


(5) 





512 


(2) 





136 


(2) 


present work 







097 


(1) 





513 


(3) 





134 


(7) 


[32] 
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Table 3 

Thermodynamic properties in the hquid phase of the binary Lennard- Jones mixture a a = o'b, 
£a = £b, ^ab = 0.75. The numbers in paranthesis denote the uncertainty in the last digit. 



T A 


Aa 


fj-B 




^B 


o* 


h* 


8%. 

HT 


V dp* 


1 


1, p^* = 0.04 
















0.05 


-4.73 (1) 


-3.87 (1) 


2.0 (1) 


1.3 (1) 


0.6892 (6) 


-5.637 (7) 


0.29 (1) 


0.13 (6) 


0.10 


-4.302 (8) 


-3.908 (8) 


1.8 (1) 


1.2 (1) 


0.6744 (5) 


-5.428 (6) 


0.31 (2) 


0.14 (7) 


0.15 


-4.140 (8) 


-3.931 (8) 


1.9 (1) 


1.3 (1) 


0.6618 (6) 


-5.262 (6) 


0.37 (2) 


-0.02 (7) 


0.20 


-4.042 (8) 


-3.940 (7) 


2.1 (2) 


1.5 (1) 


0.6506 (7) 


-5.11 (1) 


0.47 (2) 


-0.3 (1) 




1.15, p^* = 0.08 














0.05 


-4.641 (8) 


-3.226 (6) 


3.0 (2) 


1.6 (1) 


0.587 (1) 


-4.96 (3) 


0.84 (4) 


-1.5 (2) 


0.10 


-4.181 (5) 


-3.270 (4) 


2.7 (2) 


1.3 (1) 


0.559 (1) 


-4.68 (3) 


1.03 (6) 


-1.9 (2) 


0.15 


-3.952 (6) 


-3.294 (4) 


3.8 (3) 


1.6 (1) 


0.534 (2) 


-4.45 (5) 


2.0 (2) 


-4.6 (5) 


0.20 


-3.833 (6) 


-3.306 (3) 


4.4 (4) 


1.3 (2) 


0.497 (3) 


-4.17 (7) 


3.9 (4) 


-9 (1) 
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Table 4 

Vapour- liquid equilibria of the binary Lennard- Jones mixture a a = cb, ea 
The numbers in paranthesis denote the uncertainty in the last digit. 



= £b, iAB = 0.75. 



XA 


UA 


Pa 


n'* 
P 


n"* 
P 


h'* 


h"* 


1 = 


i 












0.05 


0.275 (2) 


0.0346 (1) 


0.6881 (6) 


0.0424 (1) 


-5.638 (7) 


-0.524 (2) 


0.10 


0.390 (2) 


0.0403 (1) 


0.6745 (5) 


0.0513 (2) 


-5.428 (6) 


-0.611 (2) 


0.15 


0.438 (2) 


0.0436 (1) 


0.6627 (6) 


0.0568 (2) 


-5.262 (6) 


-0.669 (3) 


0.20 


0.472 (2) 


0.0467 (2) 


0.6526 (8) 


0.0627 (3) 


-5.11 (1) 


-0.735 (4) 




1.15 












0.05 


0.158 (1) 


0.0722 (2) 


0.584 (1) 


0.0911 (4) 


-4.95 (3) 


-1.064 (5) 


0.10 


0.248 (2) 


0.0806 (3) 


0.559 (1) 


0.1051 (6) 


-4.68 (3) 


-1.158 (7) 


0.15 


0.308 (2) 


0.0890 (3) 


0.543 (2) 


0.124 (1) 


-4.49 (5) 


-1.33 (1) 


0.20 


0.341 (2) 


0.0958 (4) 


0.528 (5) 


0.144 (1) 


-4.32 (7) 


-1.50 (1) 
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Table 5 

Vapour-liquid equilibria of the binary Lennard- Jones mixture ga = £a = ^B^ ^AB = 0.75 
for the state point T* = 1.15 and xa = 0.05 from independent simulation runs. The numbers in 
paranthesis denote the uncertainty in the last digit. 



VA 


P% 




p"* 


h'* 


h"* 


0.160 (1) 


0.0732 (3) 


0.588 (1) 


0.0931 (5) 


-4.97 (3) 


-1.081 (7) 


0.159 (1) 


0.0727 (2) 


0.584 (1) 


0.0921 (5) 


-4.96 (3) 


-1.071 (7) 


0.159 (2) 


0.0720 (3) 


0.583 (1) 


0.0910 (6) 


-4.94 (3) 


-1.060 (8) 


0.158 (1) 


0.0722 (2) 


0.581 (1) 


0.0911 (1) 


-4.95 (3) 


-l.OGl (5) 
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Table 6. Thermodynamic properties in the hquid phase of the ternary Lennard-Jones mixture = ctb = ac, £b = O.TSe^, 
ec = O.lbsA, ^AB = ^AC = ^BC = 1 at the temperature T* = 1 and the pressure p'*=0.2. The numbers in paranthesis denote the 
uncertainty in the last digit. 



Xa 


Xb 


jj^A 


i^B 




V*A 


V*B 


vc 


P* 


h* 




(dh*\ 


0.972 





-3.61 (1) 




-1.69 (1) 


1.2 (1) 




2.3 (1) 


0.7114 (4) 


-5.50 (2) 


0.202 (7) 


0.43 (4) 


0.724 


0.232 


-3.890 (8) 


-3.869 (8) 


-1.79 (1) 


1.4 (1) 


1.8 (1) 


3.1 (2) 


0.6663 (6) 


-4.79 (3) 


0.35 (2) 


0.05 (7) 


0.432 


0.492 


-4.282 (5) 


-3.183 (5) 


-1.954 (9) 


1.1 (1) 


1.8 (1) 


4.0 (2) 


0.5876 (8) 


-3.83 (6) 


0.74 (4) 


-0.7 (1) 


0.312 


0.590 


-4.505 (3) 


-3.007 (4) 


-2.085 (8) 


0.6 (1) 


1.6 (1) 


5.1 (3) 


0.533 (7) 


-3.33 (4) 


1.31 (7) 


-1.8 (2) 



Table 7 

Vapour-liquid equilibria of the ternary Lennard-Jones mixture aA = ctb = crc, £b = 0.75eA, 
ec = 0.15eA, ^AB = ^AC = Cbc = 1 at the temperature T* = 1. The numbers in paranthesis 
denote the uncertainty in the last digit. 



XA 


XB 


VA 


VB 




P'* 


p"* 


h'* 


h"* 


0.972 





0.245 (3) 





0.1978 (8) 


0.7111 (4) 


0.1893 (8) 


-5.50 (2) 


-0.297 (7) 


0.724 


0.232 


0.205 (3) 


0.168 (2) 


0.207 (1) 


0.6680 (7) 


0.211 (1) 


-4.79 (3) 


-0.510 (8) 


0.432 


0.492 


0.169 (4) 


0.369 (3) 


0.212 (1) 


0.593 (1) 


0.248 (2) 


-3.84 (6) 


-0.76 (1) 


0.312 


0.590 


0.147 (5) 


0.457 (3) 


0.200 (1) 


0.532 (1) 


0.250 (2) 


-3.33 (9) 


-1.05 (2) 
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Fig. 3. 
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Fig. 4. 
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Fig. 5. 




